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ABSTRACT 


A fast, accurate and portable data reduction system 
was developed for the XR-3 Captured Air Bubble testcraft 
being evaluated at the Naval Postgraduate School, Monterey, 
California. The system consists of four units. A magnetic 
tape recorder is used for data play back. A signal condi- 
tioner unit with a built-in analog-to-digital converter was 
developed and Poeused Stomea lear amplify, sum and further 
prepare the data for transmission to either a strip chart 
recorder or a digital X-Y plotter through a Monroe 1880 
ealculator. 

Preliminary use of curve fitting techniques are dis- 
cussed; calculator programming and the various problems 
and solutions encountered in the development of the system 


are described. 
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I. INTRODUCTION 


Associated with most experimental research projects 
is the requirement to obtain, evaluate, interpret, and 
analyze large quantities of data in order to draw meaning- 
pimeconclusions. A manual data reducticn® metho@ is boch 
time consuming and prone to error. On the other hand, an 
automatic data reduction system is more accurate and can 
Significantly shorten the man-hours required to reduce 
raw experimental data into a form which can be quickly 
utilized by the project engineers. 

An automatic data reduction system was needed for the 
XR-3 Captured Air Bubble (CAB) testcraft being evaluated 
by the Naval Postgraduate School (NPS), Monterey, Califor- 
nia. This craft is being operated at Lake San Antonio, 
approximately 100 miles from the school. Test data are 
recorded on a 14-track Pemco, Model 120B, magnetic tape 
recorder. The method used to reduce the data has con- 
Sisted of transporting the tape unit back to NPS, using 
a digital voltmeter and/or a strip chart recorder to ob- 
tain a visual data presentation, and then making hand 
merformance plots. The ratio of data acq@qmisition time to 
reduction time was about 1] to 5. The need for a faster 
System arose from the decision to acquire a facility that 
would permit the test operators to remain overnight at the 


test site. In Order to maximize the efficiency of 





Operating on-site over a several-day period, the results 
of each day's testing must be reduced and analyzed before 
a schedule of tests can be prepared for subsequent test 
runs. Therefore, the time to reduce the data must be 
Shortened significantly. 

The automatic data reduction system A enesesiyad in this 
paper fulfills these objectives. It provides for a better 
utilization of operating time at the test site by provid- 
ing a quick, accurate and portable means of evaluating the 


performance of the testcraft. 
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If. BACKGROUND 


A. NATURE OF THE PROBLEM 

Both the data acquisition system and the data reduc~ 
tron methods for the XR-3 have been constantly simproved 
Since the testcraft arrived at the Postgraduate School in 
March, 1970. Initially, performance data were limited to 
thrust, velocity and plenum cavity pressure. These para- 
meters were hand recorded and hand reduced to produce 
mough Studies of thrust versus velocity and optimumeplenum: 
pressure. However, in March, 1973, a new data acquisition 
system was installed [Ref. 1] in order to test and evaluate 
the performance of the XR-3 with respect to various other 
Operating parameters such as bow and stern seal pressure, 
turning speeds, turning rates and roll damping. This is 
mie present system and the one for which the automatic 
data reduction system was to be designed and developed. A 
brief description of the data acquisition system can be 
found in Appendix A. 

The purpose of the new automatic data reduction system 
was to provide a quick and efficient means of transforming 
maw tape-recorded performance data from tests of the XR=-3 
into a finished graphical presentation capable of being 
used for evaluating the effects of varying operating para- 
INeters. The system was to have the capability of producing 


a printed record of the data converted into physical 


eal 





engineering units, such as knots for velocity vice the 
millivolts (mV) output from the tape recorder. The system 
had to be portable, in order to be moved in and out of the 
data reduction facility, and also small and capable of be- 
ing operated ina limited space. The equipment had to be 
Simple and easy to use, requiring few inputs from the 
Operator. One of the major considerations was that the 
System had to be expandable, which would allow for the 
implementation of future requirements. 

The presentation of the results was considered an 
mmmoortant factor. Three-digit data acctiracy was sufficient 
both for plotting and for the printed output since the 
raw data are accurate only to this order of magnitude. The 
data would require being scaled twice, once for the plot, 
changing millivolts to inches on the graph, and once for 
preanting, if desired, changing millivolts to physical unites. 
The system also had to be capable of performing operations 
such as summing inputs from two or more channels of the 
tape recorder. One such application is in determing the 
total drag on the testcraft at a constant speed. In this 
Bendition the drag equals the thrust, which in turn equals 
Ene sum of the thrust from both the port and the starboard 
engines. 

All of the phenomena being tested occur at a fairly 
How Lrequency, no greater than about 0.5 Hz. For this 


reason the sampling rate need not be much greater than 


1b 
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one HZ according to Shannon's sampling theorem (also Known 
as Kirchhoff's sampling theorem). However, to allow for 


future expansion a faster sampling rate would be needed. 


B. MEETING THE OBJECTIVES 

A Monroe model 1880 scientific programmable printing 
calculator and a Monroe model PL4 digital plotter were 
chosen to be the basic components of the new data reduction 
system (see Figure 1). Because of its small size and rela-~ 
tively large instruction set, the calculator seemed best 
Suited to perform any manipulations necessary to prepare 
the data for the plotter and also to provide a paper tape 
printout of the data converted to the proper physical units. 

An analog-to~-digital (A/D) converter was necessary to 
transform the analog tape recorder Set to a digital form 
Compatible to calculator input. A contract was let to 
mamta Cruz Engineering of Santa Cruz, Calitornia in April, 
1973 to design and build an interface which would make it 
possible for the calculator to select a tape recorder 
channel, then convert the signal from that channel to a 
digital input to the calculator. Calculator software (pro- 
gram) would convert this input into the various forms 
necessary to accomplish the remaining tasks. A delivery 
date of June 1973 was set; however, an extension was re-~ 
quired until late July 1973. 

As an intermediate method of data reduction while 


waiting for the interface, various methods of curve fitting 


es 





Poor Od pie DU e 





ZJOFETNOTED OSs 


SIERRA ing Me 


“OR oeceeceeuagcs. 
eens 


ES mA ROE boty 
UKAMIAS er 
Bes a. oe Hnerra 

OD i np 

RAAB se, Kies, 


VOR uO Conn 


j in WS ; a ; . a 4 ; 
TT Pe ree 


"T sanbtg 








14 





and plotting were developed on both the calculator and on 
the IBM 360/67 digital computer available at NPS. These 
methods are discussed in Section III. Because the calcu- 
lator solutions were temporary in nature and presented 
Many problems, the solutions to some of which were not 


cbhtained, only the significant portions will be discussed. 


i 





III, CURVE. FE Liitye 


A. GENERAL METHODS 

There are various methods available to engineers for 
fitting data to an algebraic relation that lends itself 
to easy manipulation and plotting and that can have a 
physical significance. The methods explored for applica- 
tion to the data reduction system were the following: 

1. Piecewise - polynomial interpolation method 

2. Difference method 


3. Least-squares approximation by orthogonal poly- 
nomials method 


4. Least-squares regression method 
1. Piecewise-Polynomial Interpolation Method 
The cubic spline method of piecewise-polynomial 

interpolation was not satisfactory for use by either the 
calculator or the IBM 360 computer for two reasons. 
First, the method required knowledge of the slope of the 
eurve at the initial and final points. Phis information 
is not known and only a gross estimate could be achieved. 
Secondly, a spline fit requires the curve to pass through 
each of the data points. Recording error and signal 
noise are not smoothed by this method and as a result an 
inaccurate model of the physical event results. A spline 
fit is best suited to determine polynomial approximations 


for exact data whose slope has discontinuities [{Ref. 2]. 
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2. 


Difference Method 


The difference method likewise does Not Lender t— 


self to computer application; however, it is an excellent 


tool for determining the degree of fit required for a 


least-squares type method. When the difference. method is 


employed, 


the data must be gathered such that the depen- 


dent variable is obtained at equal increments of the in- 


dependent variable 


A difference tabke is formed 


taking iterations until an essentially comstant difference 


is reached throughout a given iteration. The order of the 


polynomial to give the 
of iterations to reach 
That is, 


then the model will be 


if the second 


best fit is the same as the number 
this constant difference [Ret oie 


iterations produce a constant value, 


2 
= + = 2 
Of “the form. yv ay a,x a,x 


The coefficients of the polynomials are solved for by 


uSing 


where 


Yo 


the Gregory-Newton interpolation formula: 


+ x (Ay) + 


r(r-1) (A*y)) r(x-1) (r-2) (A*y,) 


+ 


a! on 


r(r-1),---, (ron) (Ay) 


n! 


is the increment between SucceSSive values of the 
independent variable x, and 


f 


are the initial values of the variables 


IL 





EXAMPLE 1 DIFFERENCE METHOD 

As an example of this method, the data for a third 
order equation y= 5 + 2x? are presented in Table I. 
Since the constant difference occurs in the third itera- 
tion, the desired fit will be given by a third order 
polynomial. In addition, since the third i enero re- 
sulted ina difference value with zero variance, the fit 


will be exact. 
TABLE I 


Successive Differences for Difference Method of Curve 
meeting. 


Successive Dwtrtecrences 





lst 2nd Se 4th 
Given Data Iteration Iteration| Iteration Iteration 
x y Ay a4 


4 
0 sMAyES 16 | A*y = Sel A comm aT a 
0 , 6 =e ee 
2 Ze Ay,= 12 A Y,= 1924 Y,7 96 A y.= re) 
4 133 | Ay,= 304 A“y = 288 Amy, 96 
6 437 | Ay,= 592 A*y4= 384 
B 1029 | Ay,= 976 
10 2005 


a 
2 2! 


M< 
It 
) 
ay 
NO 
Ma 
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which is the exact equation from which the data were de- 
rived. To determine the order of a model from data which 
is not exact, each column of the difference table must be 
inspected for an essentially equal distribution of plus 
and minus signs. The number of iterations to that point 
determines the order of the resulting polynomial fit. 

The remaining procedure is as discussed in the preceding 


example. 


3. Least-Squares Approximation sy @menoganat 


Polynomials 

The. method of least-squares using orthogonal poly- 
nomials could not be used as a calculator solution to the 
problem since more storage is required by the method than 
is available in the calculator. However, it does provide 
a very efficient and accurate solution when used on the 
IBM 360/67 computer. 

The method consists of solving for successively 
higher order orthogonal polynomials, such that the sum of 


which satisfies the requirement that 


N 
2 
> bi(x ) =] pix 8 where 
n n 

n=l 

f(x) = exact solution 

p(x) = polynomial approximation 

N = order of the polynomial 


is as small as possible [Ref. 2]. A computer program was 
developed for use on the IBM 360/67 in which the X and y 


data points, the number of points and the degree were 


ED 
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entered. The output consisted of two parts: The first 
contained: 
1) The coefficients of the fitting polynomial 
2) An estimate of the error of each coefficient 
3) Sum of the square of the deviations 
4) The F-ratio 
5) The Chi Square goodness of fit Statistic 
6) Degrees of Freedom of the data, and 
7) A list of the data points, the evaluated poly- 
nomial at each value of X, and the differences 
between the two y values. 
Batt two 1S a plot of: 
1) The data points, 
2} the fitted curve, and 
3) the normalized error, 
The computer program and sample output can be found in 
Appendix B. 
4. tLeast-Squares Regression Method 
The least-Squares regression method of data fitting 
is the most familiar and the most widely used by engineers. 
It is used by hypothesizing the existence of a polynomial 
that will best fit the data as in the case of orthogonal 
polynomials. The degree of the fitting polynomial must be 
determined by another method or by analysis of the results 
moeorder to determine the "best" fit. In the case of an 


Nth order fit, the fitting polynomial has the form: 


2 3 N 
Va oe er Ce ee eee Tk 


20 





eit 





Additional equations are needed to solve for the coeffi- 


client a, b 


plying the equation by 


variable on both sides 


+b ), x 
BY Se 155 ye oo 
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x yY 
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can be written: 


. (2) 


(3) 


EXAMPLE 2 Least-Squares regression method 
As an example, the data of example one will be fitted 
to a third order polynomial. The data in Table I can be 


used to generate the data in Table II. 


TABLE II 


Least-Squares Data 


x2 oe xy a xy 2 6 

x 

5 a 
| 16g 64 


Sa eae ee 4096 
é 16656 
g|1o20|_64| 8232| 512 | 65856 | 4096 | 52684e |32768| 262144 
Lo 100000 


E [3013630 220 |31478| 18600 (2843500 (1566491263 4920 (S160 0G)| Zz ce 


Substituting the data from Table II into the Equation (1) 


3630 = are! 30b + 220c + 1800d 

31478 = 30a + 220b + %1800c + 15664a 
284300 = 220a + 1800b + 15664c + 516004 = 
2634920 = 1800a + 15664b + 51600c + 4129604 


Solving (4) for the coefficients gives 
a= 3,b = 0, ¢ = Ofrd = 2; 


therefore the resulting equation is 
3 
y = Set 2x 


which is the exact solution to the data being fitted. 
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B. CALCULATOR APPLICATIONS 

Experimentation with representative data on the IBM 
360/67 computer showed that a polynomial of seventh order 
was required to best-fit a curve smoothly through the 
data points. This was also verified using the difference 
method. The least-squares regression method of curve fit- 
ting was implemented on the calculator solving for a 
seventh order polynomial. Various problems were en- 
countered while trying to solve Equation (1) for the un- 
known coefficients (a,b,c, ..., n) because the A matrix 
was highly ill-conditioned. Two methods were used to try 
to solve the equation. The first was a direct approach 
uSing triangular factorization with forward and back sub- 
Stitution. The second was an iterative method using a 
Gauss-Seidel iterative scheme with and without successive 
Over-relaxation (SOR). Examples of these methods can be 
found in reference 4. 

The triangular factorization did not produce accept- 
able results because of the loss of significant digits. 
The calculator solution differed so extremely from the 
computer solution and the plot of the resulting polynomial 
equation using the calculator solution bore so little 
resemblance to actual performance plot that other methods 
Naa to be tried. Scaling of the input data had no correct- 
ing effect on the results using this procedure. 

The Gauss-Seidel iterative method occaSionally would 


‘ 


provide the proper coefficients for the fitting polynomial; 
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however, this was also an unsatisfactory method as a 
general solution to the problem. See Figure 2 for the 
flow chart of the calculator program implementing the 
Gauss-Seidel iterative method. The initial guess that 
was used to arrive at the proper solution was extremely 
critical. An improper guess would lead to instability 
and the coefficients would diverge from the proper solu- 
tion. Even when stable and converging, rie time required 
for the calculator to reach a solution was excessive. 
An iteration required about 20 seconds and, depending on 
the initial guess, a problem could require over 200 itera- 
tions to arrive at a solution which would closely approxi- 
Mate the physical phenomenon. The introduction of a 
relaxation factor helped; however, like the initial guess, 
the choice was critical and no way was developed to aid 
in the choice of either an initial guess or a relaxation 
merectOr , 

The time problem and the initial guess and relaxation 
factor problems caused this method to be unacceptable as 
a general solution to the problem of curve fitting on the 


Monroe 1880 calculator. 
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depress JUMP, O, O 
load cards 1-1 and 1-2 
place SENSE SWITCH up 
depress JUMP, SYM, 1, RESUME 


SYM 1 enter X, depress RESUME 
enter Y, depress RESUME 
NO 
SYM 2 enter initial guess for Gauss Seidel iteration 
depress RESUME 


SYM 3 











depress JUMPjm=OoF 0 
load cards 2-1 and 2-2 

place SENSE SWITCH up (to iterate) 
depress JUMP, SYM, 6, RESUME 


SYM 6 enter relaxation factor 
depress RESUME 














SYM 4 calculator iterates and prints coefficients 
a thru a_ of a+ a.xt+ a.x7+...4+ ax 
every fifth iteration : d 

SYM prints present iteration then stops 

CHANGE 

SIGN 





Figure 2. Flow Chart of Calculator Curve FPCEINng Rotttane. 
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check 
present 
solution? 






NO 











depress JUMP, SYM, 9, RESUME 


enter A, depress RESUME 
calculated y is printed 









continue 
iteration with 
present values? 







depress JUMP 
SYM, 4, RESUME 







to start with 
new initial 
guess? 


depress JUMP, 
SYM, © RESUME 
emter new 
guess 

depress RESUME 











NO 





TO 
PLOT 
ROUTINE 


Pagunes . (continued) 
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IV. DATA REDUCTION SYSTEM 


A. SYSTEM COMPONENTS 

The data reduction system is composed of the following 
five major components: 

1. Signal selector and conditioner unit 

2. A/D converter and calculator interface module 

3. Monroe 1880 calculator 

4. Monroe PL4 digital plotter 

5. Hewlett Packard model 7100B strip chart recorder 

The data in the form of an analog signal are played 
back from all 14 channels of the Pemco magnetic tape re- 
eerder. These signals go into the signal selector and 
conditioner unit from which the operator chooses the path 
that he wants each of the signals to follow. He can send 
up to eight signals to the A/D converter module also 
located in the signal selector and conditioner unit, which 
in turn feeds digital raw data to the Monroe 1880 calcula- 
tor for manipulation and plotting, or he can send the raw 
analog data to the multi-channel strip chart recorder. 
This flexibility is essential to the operation of an 
efficient data reduction system. 

ion olondl  Seleceor andsecond? tiene tee minr 

The signal selector and conditioner module was de- 


Signed to allow for operator ease in selecting, monitoring, 


eu 
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and routing signals from the tape recorder to the other 
components of the data reduction system. All channels of 
the recorder feed into the module through the first 14 
terminals of a CO-ORD Switch 10 x 15 pin, 2-deck Matrix 
Program Board. By selecting the desired output channel, 
the operator can filter, amplify and/or sum signals uSing 
any of nine filter/amplifier circuits and a summing cir- 
cuit on the 10 outputs of the matrix (see Figure 3). The 
Signal voltage can be monitored at any of 24 points uSing 
a built-in Dee Sil model DM 2000AR Digital Panel Meter and 
a 24 position rotary selector switch. The monitoring 
points include each of the inputs from the 14 tape re- 
corder channels, each of the outputs from the nine filter/ 
mmelifier circuits and the output from the summing circuit: 

Signal conditioning was needed to filter out the 
extraneous high frequency noise from the tape recorder 
Signals which originated from the XR-3 onboard sensors. 
The noise resultS principally from vibrations of the test- 
craft during normal operations; consequently each of the 
Signals has about the same characteristic frequency noise 
as shown in Figure 4, The filters used are a Butterworth 
low pass type [Ref. 5] with a cutoff frequency of 0.5 Hz. 
The electronic circuit diagram and the signal to frequency 
response curve are shown in Figure 5. The circuit was 
constructed uSing a Philbrick Nexus SQ-10a amplifier. 

The output from the filter is then amplified using 


a Burr Brown 34403 amplifier circuit. The one-volt 
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meagure 3. Signal Selector and Conditioner Unit. 


29 








7 ~~ a . 
' > —- 
=> ¢€ 
‘ 
>< «€ 
; a a <a 
o- 





x 
Port Thrust 


Starboard Thrust 


Roll 


Pitch 


Yaw 


Yaw Rate 


melocity 


Rudder Position 


Forward 
Velocity 





LO 20 30 40 50 60 70 80 90 


a Time (m seconds) 
1 ecm = 0.1 volt 


Figure 4. Sensor Output Signals. 
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Maximum signal from the tape recorder must be amplified 
to a five-volt maximum signal for use with the A/D con- 
verter. This is to ensure the best possible data 
accuracy when uSing the A/D converter and calculator in- 
terface module. Since each of the sensor/amplifier/ 
Signal conditioner circuits on the testcraft does not 
Maintain a null output for a zero signal input, a zero 
offset voltage circuit was included. Inasmuch as a full 
scale input does not always have an exactly 1.0 volt out- 
put, a ee gain resistor was placed in the feedback 
circuit of the amplifier to allow the aperator to adjust 
the gain to have a 5.0 volt output for a full scale input. 
See Figure 6 for the amplifier circuit. 

Two of the filter/famplifier circuits are added 
then halved through a summing circuit uSing another Burr 
Brown 3440J amplifier (see Figure 7). This circuit is 
used in obtaining total thrust from the individual port 
and starboard thrust components. 

Electrical power is supplied to the Burr Brown and 
Philbrick amplifiers from a, Philbrick/Nexus NPS-300 +1l15v 
-15v power supply. Zero offset voltage 1s Supplied by a 
SRC MOD 3564 power supply and the Datel Digital Panel 
Meter is powered by a Acopian MOD 10375 10v dc power supply 
reduced to 5v dc by a National Semiconductor LM340K-7805 


5v dc regulator. 
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ie Analog-to Opie Uc GU Converter and Calculator 


Interface Module 

The analog-to-digital converter and calculator in- 
terface module, also called the A/D interface, was pur- 
chased from Santa Cruz Engineering in late July 1973. Many 
modifications were required to permit it to be used with 
the Monroe calculator and the other components of the data 
reduction system. These modifications are discussed in 
Section V. ‘ 

The A/D interface is selected through connections 
on the front of the signal selector unit. AS Many as eight 
Signals can be connected and generally the outputs from the 
filter/famplifiers are used as the inputs. The calculator 
program selects which input signal is to be converted by 
sending a three-digit octal code from the calculator IX 
register to a TTL/MST 74164, S-bit, serial to parallel con- 
verter in the interface. The first three bits are then 
sent to a RCA COS/MOS analog multiplexer CD4051A for the 
channel selection. The multiplexer receives the input 
Signals from all eight input channels. It acts as a binary 
controlled analog switch by selecting which channel is to 
be outputted to the A/D converter. The A/D converter 1S 
then triggered by a TTL/MONOSTABLE 9602 multivibrator. 

The A/D converter is a high speed, highly accurate 
and stable Analogic MP 2410 unit housed in a small 4 inch 


x 2 inch x 0.4 inch package. It is capable of ten-bit 
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resolution; however, only eight bits are used for com- 
patibility with the Monroe 1880. It is capable of ranges 
of -10V to +10V, -5V to +5V and OV to +5V with output in 
one's complement, two's complement, offset binary or 
unipolar binary. The OV to +5V range using unipolar 
binary was selected for the data reduction system appli- 
cation. The output from the A/D converter then goes into 
a TTL/MSI 7491, eight-bit, serial-in serial-out, shift 
register. Then upon another command from the calculator, 
the binary Signal is transferred from the shift register 
in the interface to the IX register of the calculator for 
further processing. 

Timing for the interface is provided by the A/D 
momnverter for the internal function and by the calculator 
for data transfer. Power is provided by an Analogic 
AN3001 modular power supply outputting +15V, -15V and 
mo aC. Timing and transfer logic is provided by FAIR- 
CHILD TTL/LSI 7400 Quad two-input NAND GATE and 7410 
triple three-input NAND GATE integrated circuits. The 
logic diagram for the A/D interface can be found in 
Figure 8. 

Semonroe 1880 Ccalculator 

The Monroe Model 1880 Scientific Programmable 
Printing Calculator was chosen as the heart of the data 
reduction system because it combines the features of a 


simple-to-operate but powerful scientific calculator with 
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an extensive problem-solving capability. The model being 
used provides 1024 memory locations for program steps (in- 
Bt mileemons 128 main data registers and 10 temporary 
(scratch pad) registers (both for data storage). The main 
data registers can be programmed to hold an additional 
1024 instruction steps instead of data if a larger program 
is necessary. Data is stored to an accuracy of 13 signi- 
Picante digits. 

The operations of the calculator fall under six 
general categories: Inpiae: sepa Storage, control, 
execution and arithmetic computation. All of these opera- 
tions can be programmed from the keyboard either through 
individual function keys or through an ENTER CODE key 
which allows the three-digit octal code of the instruction 
to be entered into memory (See Figure 9). 

A program can also be stored using a magnetic card 
device which is built into the calculator and enables pro- 
grams to be both read into the calculator and written on 
to small magnetic cards each capable of holding 512 program 
steps. Data from the main data registers can also be 
written on the cards and conversely data from the cards 
can be read into the main data registers. Data may be in- 
putted to and outputted from peripheral devices through 
the IX register and an external 36=-pin connector to the 
peripheral. Programs can be printed on paper tape with 
the listing containing the instruction number, the three- 


digit octal code and a symbol to identify the instruction 


37 





1 we OG Ae i oe 7 


OOTUOW 


wo 


gS aM 


SoG O16 Fal 





AS. Wy wwweFoFEES 


38 





if one is applicable. Program steps can also be inserted, 
deleted or changed without having to re-enter the entire 
code. This feature is a tremendous aid in correcting 
errors in a program. 

Control functions include such items as sitting 
flags, printing data, and setting the position of the 
decimal point. Program execution is extremely easy with 
the Monroe 1880. Symbolic program addressing makes it 
possible to locate a program anywhere in program memory. 
It permits branching or jumping to any sips in a program 
or to the ey ere of a subroutine without having to 
know the absolute storage location address. Branching re- 
turns to the point of the branch after the execution of a 
reserve; a jump does not. Decision making is made simple 
with the use of a sense switch, a programmable and a key- 
board flag and through the test of a positive, negative 
or zero value of the contents of the main entry register. 
Tests of these devices are used to control branches or 
Jumps to other parts of the loaded program. 

The arithmetic functions performed by the calcula- 
tor include such standard keyboard items as add, subtract, 
multiply, divide, square root and reciprocal. in add1t2on- 
other keyboard items such as sine/cosine, are sine/arc 
cosine, logarithm, statistical summation and factorial are 
included. Many more functions and controls are available 
through enterable ingtruction codes for which keys are not 


available. These include most of the codes for transferring 
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information to peripheral devices, program control, regis- 
ter manipulation and other arithmetic functions such as 
absolute value and arc tangent. A complete set of instruc- 
tion codes and operating instructions can be found in 
Refs. 6, 7, and 8. 
4. Monroe PL4 Digital Plotter 

The Monroe PL4 Digital Plotter is a peripheral de- 
vice used in association with the Monroe 1800 series calcu- 
lators. Its purpose is to display graphically data which 
is passed to it as a result seeeee eee execution. The graph 
origin may be set anywhere on the plottine surface through 
the use of a calculator scaling routine which also formats 
the data for transfer to the plotter. The plot subroutine 
discussed in Section IVB and shown in Appendix D, determines 
whether the plotter is “busy” (actually plotting) or is 
ready to receive data. If busy, the program loops until 
the plotter is ready at which time it passes the pen command 
and the data for plotting. 

Instructions and data are passed to the plotter 
through the calculator IX register. The four commands 
which can be issued by the calculator are presented in 


Table III. 


TABLE IIT 


Command Codes for Monroe PL4 Plotter 


CODE COMMAND 
160 X data to be passed 
iis al Y data to be passed 
OZ PEN UL 
163 PEN DOWN 
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The plotter responds to the PEN UP and PEN DOWN commands 
only when the REMOTE switch is depressed (see Figure 10). 
In the manual mode of operation, when the REMOTE switch 
is not depressed, the pen responds to the commands issued 
by the operator. These commands are: 

WRITE (pen down) 

EEN UP 

MOMENTARY (spring loaded up, pen down while depressed) 
and are issued through switches on the plotter control 
panel. After the X data or the Y data instruction has been 
passed to the plotter, the calculator transfers the data 


from memory to the IX register for transmission to the 


plotter. See Table IV for further plotter characteristics. 


TABLE IV 


Monroe PL4 Digital Plotter Characteristics 


Plotting Area 10 inches x 10 inches 
Ballpoint Pen Width O07 COS einen 
Resolution (X and Y axis) PrOoOGwecounts 
Repeatability O.l percent of -full scale 
Accuracy 
Static O.2 percent of full scale 
Dynamic 3 inch retrace will overlay 


trace within 0.035 inches 
between center lines 
Dimensions 


Width 17 inches 
Length 21 inches 
Height 767 5S-sinehes 
Weight 48 pounds 


aA Hewlett Packard Model 7100B Sie Oe cilongemnerc Order 


The Hewlett Packard 7I0O00B Strip Chart Recorder 
(Figure 11) has two independent pen drives and input modules 


which can accommodate analog input up to 100 V. The chart 
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transport system has 12 selectable speeds and uses a 
Standard chart roll with a 10-inch writing width. The 
pens move horizontally across the top of the chart each 
capable of a full 10-inch travel. Specifications are 


listed in Table V. 


TABLE V 
HP 7100B Strip Chart Recorder Specifications 


Response time 0.5 seconds full scale 
Chart Speeds Led iny he0 .1 ie 2 . Oy 
245 ay mn 
O. 1, O27 ORer 1 , 2 snes 


Input Range Variable from lmV to 100 V 
Input Resistance 1 megohm 
Dimensions 

Width 16-3/4 in 

Height 8-11/16 in 

Depth 7-1/4 in 

Weight 20 pounds 


Bis CALCULATOR PROGRAMMING AND DATA TRANSFER 
1. Memory 

The calculator memory is divided into four pages. 
Each page is then further divided into columns with each 
containing 32 registers, numbered in octal from 00-37. 
The first page, Page 0, is user's programmable memory and 
Pencains four columns. Each,column can hold 256 instruc- 
tion codes giving a total capacity of 1024 instruction 
Storage locations. The second page (Page 1) is data 
register memory also containing four columns for a total 
of 128 data registers used for main data storage. The 


third and fourth pages contain 16 columns each consisting 
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of working storage and read-only memory (ROM). The ROM 
contains the permanently stored programs that perform the 
keyboard functions. The last two columns of the fourth 
page (Page 3) containing the working storage, and all of 
the first two pages are made up of random access memory 
(RAM) which can be written into and read from. Informa- 
tion stored in Pages 0, 2 and 3 is treated aS a program 
instruction. A single instruction (or program step) 
occupies one byte (8 bits). Eight Single instructions can 
be stored one item per register. The format is shown in 
Figure 12. Each byte contains two Binary Coded Decimal 
(BCD) digits except the Sign position. 

The calculator contains three programmable control 
registers, the Entry (E) Register, the data Address Regis- 
ter (DC) and the Index (IX) Register. The entry register 
is the primary data register. Inputs from the keyboard are 
put in the E register. Output data are formatted there, 
then transferred to a print buffer and at the end of any 
calculation the E register contains the result. The entry 
register is formatted the same as a data register with the 
addition of an overflow and an underflow digit. The address 
(DC) register contains 14 bits and contains the address of 
a specific location in memory. It is used with various 
Macro instructions for defining the address to be used for 
data or code manipulation. 

The index (IX) register is a general purpose 8 bit 


register. It is divided into two 4 bit parts, IH (index 
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high) and IL (index low). The IX is used as a buffer 
Begister for the input and output (170) of S@nEormation to 
peripherals through the I/O bus. All information trans- 
fer between peripherals and the calculator occurs through 
the IX register. 
2. I/O Lines 

The data output line from the calculator to the 
moO bus is BYOT (byte, output, pin 25 on the connector). 
Data from the IX register is transmitted on this line on 
the falling edges of the clock pulses from T8-T15 time 
(see Figure 13) during an OPCB (output control byte, code 
314) command and during an OPDB (output data byte, code 
315) command. The OPCB command is used to transmit selec- 
tion and control data of peripheral devices while the 
OPDB command will send numeric data over the output line. 

‘The data input line into the calculator from peri- 
pherals is BYIN (byte, input, pin 26). Data is.transmitted 
from the peripheral when an IPDB (input data byte, code 
305) command is executed. 

The calculator transmits two signals to the, peri-— 
Pheral device to allow the device to determine which type 


Se,command, control or data I/O, the calculator is execut- 





ing. The signal OTCT (output control time, pine 23) re- 
Mains at a logical O from T8-T15 time and DTEN (data 
transfer enable, pin 24) goes to a logical 1 when the 
calculator executes an OPCB (output control byte). If an 


OPDB (output data byte) or an IPDB (input data byte) 
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Figure 13. Data Transfer Signals. 


command is executed, DTEN will remain at a teogicaly Ostrom 
T8-T15 time and OTCT will go to logical l. 
3. A/D Interface Data Transfer 

The operator determines which of the input signals 
to the A/D interface are converted and sent to the calcu- 
lator by entering the. proper code (shown in Table VI) into 
the computer program at decimal step numbers 130 and 135 
(see Appendix D). During execution, the calculator loads 
the IX register with the code for the channel of the 
dependent (y) variable (decimal steps 129 and 130) then 
outputs a control byte (OPCB code 314) to set the multi- 
plexer to the proper channel then convert the signal. The 
next step (IPDB code 305) transfers the converted signal 
to the IX register. These data are then stored and the 
code for the channel of the independent (X) variable is 
loaded into the IX register and the process is repeated. 
After the two signals are in the calculator, they are con- 
verted from binary to BCD, scaled, and printed if desired, 


then passed to the X-Y plotter. 


TABLE VI 


A/D Interface Channel Codes 





CHANNEL CODE 
iL 340 
2 140 
3 240 
4 040 
5 300 
6 100 
7 200 
8 000 
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4. Digital Plotter Data Transfer 

The digital plotter has a slightly more compli- 
cated data transfer scheme because of itS slower response 
time. The first 42 instructions (000-041) comprise the 
Pemotter subroutine and are used to control data input to 
this device. The main program branches imto the sub- 
routine at the decimal steps 000, 004, OO8 and 012. When 
a value is to be passed, rather than a pen command, the 
value is stored in the E register and is register 27 of 
the last column. For example, when transfer of Y data to 
the plotter is desired, the program branches to step 004. 
The program in Appendix E can be followed. 

The code for Y data (161) is loaded into the IX 
register than the subroutine jumps to octal step 16. The 
contents of the IX register is then copied into another 
register. The IX retains a copy, then the calculator out= 
puts a control byte (OPCB). The program then checks the 
lower 4-bits of the IX register (IL) to be not equal to 0. 
Since the IL register contains a l, a jump to octal 
address 25 (decimal address 21) is performed. A RINS 
(read input status) 1S executed and the plotter returned 
alto the H8 (left most) bit of IX 1£ busy and the 
calculator repeats until not busy or the plotter returns a 
O. The IX register is then re-loaded with the initial Y 
code (161) destroyed by the last step then tested for a 
pen command. This test fails so the sign digit of the E 


register is passed to IL then transferred to the plotter 
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(OPDB). The address register (DC) is loaded with the 
address of register 27 of the last column by the instruc- 
tions 273 and 277. This register contains the number in 

the E register; however, data transfer is easier from 

this location. IX is loaded with the contents of the 
address pointed to by DC which is the two exponent digits. 
These are transferred to the IX register then DC is de- 
cremented by one step (two digits) and the transfer is 
repeated. The plotter receives the exponent, the four 
most Significant digits of Sere and the sign of both. 


The plotter pen then moves to the point. 
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V. A/D INTERFACE PROBLEMS AND SOLUTIONS 


Initially the A/D interface had four features which 
inhibited it from being used with the rest of the data 
reduction system. The two major problems were (1) the 
logic for data transfer and (2) the loss of the most sig- 
nificant bit from the A/D converter. The two less serious 
and more easily corrected faults were (3) the voltage eon 
version range of the A/D converter and (4} the conversion 


of binary to decimal numbers. 


A. LOGIC FOR DATA TRANSFER 


The initial design of the A/D interface assumed both 








DTEN and OCTC to remain logical O's throughout the entire 
normal time cycle instead of both going high from T8-T15 
time. It also assumed that DTEN would goyvtigh for data 
transfer and OCTC would go high for control output instead 
of missing the one time pulse train. Consequently the A/D 
interface could not interpret the calculator requests. In 
Seeger to correct this design flow, the circuit labeled "A" 
in Figure 8 was incorporated. The two signals, OcTC and 
DTEN, are compared with an exclusive OR. If the two are 
different, the output of the exclusive OR is high. The 
original OCTC and DTEN are then "AND"ed to determine if 
either signal has been missed indicating the command has 
been executed. This addition changes the OCTC and DTEN 


Signals to that expected by the rest of the interface. 


3) 7% 








B. LOS a 1 OF MOST eSIGNiEI CANT BET 

After correcting the first problem, analog signals 
could be converted to binary and transferred to the cal- 
culator. However, only seven of the eight bits were 
being passed to the IX register. It was determined that 
the most significant bit was not transferring into the 
interface 8-bit shift register. It was left "floating" 
at the output of the A/D converter then lost when the 
first transfer pulse was being issued by the calculator. 
The circuit labeled "B" in Figure 8 was added to remove 
the first transfer pulse during an IPDB and allows the 
entire 8 bits to be transferred into the IX register. 
It uses a TTL/MST 7493 4-Bit Binary Counter and four "NOR" 
gates configured to remove the T9 time pulse from the 
clock to the 8-bit shift register during the transfer of 


data from the A/D converter. 


CG. VOLTAGE CONVERSION RANGE 

The Analog-to-Digital converter was wired at delivery 
to convert signals ina range of -5V to +4.96. Since all 
the signals from the XR-3 are ina range of OV to #+l1V, 
this larger range reduced the significance which could be 
achieved by the converter. Using eight bits coded binary, 
the full scale range of the converter can be broken down 
into ~~ Or 256 parts. The smallest range capable of the 
A/D converter is OV to +4.98V so in order to realize the 


greatest possible accuracy two changes had to be made. 
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First the A/D converter range was changed to OV to +4.98V 
and second the signal condition and amplifier unit was 
developed and built to enable the tape recorder signals 
to be amplified to a 5V maximum signal. By making these 
two changes the full capability of the A/D converter is 


utilized. 


D. BINARY TO DECIMAL CONVERSION 

The calculator data format is binary coded decimal 
(BCD) characters. Four bits are used to represent a deci- 
mal digit and a number is made up of by combining BCD 
characters. The A/D interface transfers a number to the 
IX register in binary representation. In order to be used 
by the calculator, this number has to be converted to BCD 
format. This is done by calculator software by a series 
of jump instruction beginning at decimal step 48 of the 
data reduction program. IX register bit Ll (right most) 
is the most significant bit. Lf this bit is high 2s eas 
feed tO scratch pad register 0. If bit L2 is on, 64 is 
added, and so on such that if all bits are on, scratch pad 
register O totals 255. After checking all bits of the IX 
register, the scratch pad register is divided by 256. The 
result of this is to leave in the register the decimal 
equivalent Soe Ov to 1V analog signal converted by the 


Pee interface. 
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EXAMPLES OF CONVERSION: 


VOLTAGE IN A/D OUTPUT 


Oe 2lV 0100 0000 
7 O0V 1000 0000 
oe 75V 1011 0000 


OV 1 LL gt Ee a i 


CALCULATOR CONVERSION 


ye 
27/28 
(2742° 


(2742° 


55 


64/256 = 0.25 


128/256 = 0.5 
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VI. CONCLUSIONS 


The automatic data reduction system has been shown to 
provide a quick and accurate means of evaluating performance 
data from the XR-3 testcraft. It has been tested with both 
Simulated and actual tape recorder signals and proven to be 
accurate to greater than 1% of full scale. Raw data can be 
efficiently transferred into a finished graphical presenta- 
tion and a printed record of the data, converted to physical 
units, can be easily obtained. The system is superb for 
its simplicity and ease of operation. A job that used to 
require days can now be completed ina matter of hours. 

The signal selector and interface unit provides the 
capability of zeroing and calibrating the signals from the 
tape recorder. It removes the signal noise resulting from 
vibration and can pass the data to a strip chart recorder, 
Seeit can internally digitize the information for transfer 
to the calculator for further processing and/or plotting. 

The provision for expansion has been designed into 
the system both through the modification of hardware and 
software. The Signal selector and conditioner unit has 
room for internal modification and additions to satisfy 
requirements which may arise in the future, and the calcu- 
lator has over six columns of unused memory which may be 
used to increase the data reduction program to fulfill 


mreuare needs in this area. 
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APPENDIX A 


XR-3 DATA ACQUISITION SYSTEM 


The data acquisition system presently installed in the 


XR-3 testcraft is comprised of various sensors feeding 


through amplifiers and/or signal conditioners in order to 


prepare the signal for recording on a Pemco Model 120-B 


Magnetic tape recorder (see Figure 14). 


ing 


bridge transducer. 


SENSOR SYSTEM 


Sensors are presently installed to measure the follow- 


werameters ; 


On DM & WD EH 
ry 


\O 
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a0. 
a. 
are « 
a3. 
14. 
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Thrust 1S measured by a Revere USPI-150A 


Port. thrust 

Starboard thrust 
Forward seal pressure 
Aft seal pressure 
Plenum chamber pressure 
Testcraft velocity 
Water emersion height 
Pitch angle 

Pitch rate 

Roll angle 

Roll rate 

Yaw angle 

Yaw rate 

Lateral acceleration 
Longitudinal acceleration 


(30) balance 


One is installed on each engine mount 


sO as to measure the fore or aft thrust of each engine. 


forward and one aft. 


1 


Two pressure probes are installed in the plenum, one 
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Only one can be used at a time. 











Balance Bridge 
Thrust Transducer 
(Starboard Thrust) 






Balance Bridge 
Thrust Transducer 
(Port Thrust) 


Amplifier- 
Transducer 
Package 


Forward 

Seal Varian 

Pressure Differ- 
ential 
Trans- 

Aft ducers 

Seal 

Pressure 


Plenum 
Pressure 









Velocity 
Vane 
(FM Signal) 





Velocity 
Conditioner 


Pemco 
Magnetic 
Tape 
Recorder 








Height 
Sensor 


Height 
Sensor 
Conditioner 





Pitch angle 


piaphrey's Roll angle 


@-16-0101-1 Yaw angle 


Gyro Pitch rate 
Package Roll rate 


Yaw Rate 





Lateral acceleration 


Presently 
Longitudinal acceleration not used 





Figure 14. Data Acquisition System. 
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The signal from the bridge is amplified by one of ten Grant 
Model DCAB-3 amplifiers to a range between 0.0 and 1.0 volts 
eorresponding to 0.0 to 500.0 pounds thrust. Calibration 
Signals of 0.0, 0.5 and 1.0 volts are also available to 

each of the amplifiers which are found in an amplifier- 
transducer package. 

Pressures are measured by pneumatic probes connected to 
femerdan differential pressure trans@ucers jaccurate@to ur ie. 
Pressure Signals are amplified and calibrated similar to 
thrust so as to have 60 ings HERE correspond to 1.0 volt. 

Velocity is measured by a vane enclosed in a water probe 
which transmits a frequency modulated signal to a Potter 
velocity conditioner. The velocity conditioner converts 
the FM signal to a dc Signal from 0.0 to 5.0 volts corres- 
ponding to 0.0 to 40.0 knots. This signal is then reduced 
eo, 0.0 to 1.0 volt by another amplifier in the amplifier- 
transducer package. Calibration is performed by operation 
Over a measured course. 

A Humphreys Model CF18-0O101-1 gyro package measures 
pitch angle and rate, roll angle and rate, yaw angle and 
rate, and lateral and longitudinal acceleration. Presently 
all modes except acceleration are being used. Outputs are 
between 0.0 and 1.0 volt and range is + 20 degrees in pitch, 
mes degrees in roll, + 180 degrees in yaw, and + 30 degrees 


fJsec in all rates. 


D2 








Height is measured by a Western Marine Electronics 
Model LM4001 height ee with the transducer located 
eorward of the bow. Range is seliectan ieee cect or 
Bees feet, both with an accurate of lene TZ etece 
range is presently in use. The signal from the sensor 


is converted to 0.0 to 1.0 volt de by a height sensor 


conditioner. 


B. DATA RECORDER 

The conditioned signals from the various on board 
sensors are sent to and recorded on a Pemco Model 120-B 
Magnetic reel-to-reel tape recorder. The recorder has 14 
channels for electronic data plus an edge track for audio 
recording. The unit can be controlled from either the 
pilot's instrument panel or from the unit itself, located 
just aft of the pilot's cockpit. The recorder is Sor eaDme 
weighing only 100 pounds and measuring 9-1/2" x 17-3/4" 
x 26-1/2". It is easily installed and removed from the 
testcraft. 

The recorder can operate at tape speeds of 60, 30, 15, 
7-1/2, 3-3/4, and 1-7/8 inches per second and uses 
standard or precision NAB 10-1/2 inch reels and 1.0 or 


1.5 mil one-inch tape. 
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APPENDIX B 


IBM 360/67 COMPUTER DATA FITTING 


A fortran program was developed for use on the IBM 
360/67 computer to aid in calculator software design. 
The program consists of five basic parts as follows: 

Pages: initialization and contre 

Part 2: Data point plotting 


Part 3: Determination of fitting polynomial coeffi- 
cients and error criteria 


Part 4: Error normalization andgplctemng 

Part 5+: Evaluation and plotting Gletina! curve 

The plotting subroutine, DRAW, is used to provide the 
final graphical output. It is available as a callable 
subroutine and utilizes a Calcomp Plotter. 

The subroutine LSF provides the fitting polynomial 
coefficients and all the written output. It is a modifi- 
cation of the subroutine LSQOF2 and the procedures and 


algorithm for the method can be found in Ref. 2. 
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m WhO e 
e 


Om Bh WH FE 
° 


sim ° 


SYM O 


SYM l 


SYM 2 


depress JUMP, O, O 
load card "DATA REDUCTION ROUTINE" 

place sense switch up for point plot 
depress PEN UP switch on plotter 


APPENDIX C 


FLOW CHART FOR DATA REDUCTION PROGRAM 






* 
CHANNEL CODE 










a ce 
depress JUMP, l, 3 Zee we Lad 
place RUN/STEP/LOAD switch to LOAD eS Gee) 
enter code for y axis channel* Ge ecs 040 
depress LIST PROG 4 times Doe ae O00 
enter code for x axis channel* Sono oleate, 
put RUN/STEP/LOAD switch to RUN aie warecou 

Sao A Ole: 


“OMB WD EF 















depress JUMP, SYM, 0, RESUME 
set 0.00 (button left corner)on grid with ZERO controls 
depress RESUME 

set 10.00 (top right corner)on grid with SCALE controls 
enter X grid coordinate of origin, depress RESUME- 
enter Y grid coordinate of origin, depress RESUME 

pen moves to origin, depress REMOTE switch on plotter 








draw NO depress JUMP, 
axes? | SYM, 5, RESUME 


YES | (5 









enter X scale (units/0O.1 in), RESUME 
enter Y scale (units/0O.l in), RESUME 
read length of x and y axes 


depress JUMP 
SYM, 2, RESUME 


ee, 


5. jf) 





SYM 3 


SYM 5 









enter interval between graduation of X, RESUME 
2. enter interval between graduation of Y, RESUME 
axes are plotted 


1. depress JUMP, SYM, 5, RESUME 









ik enter x scale (in/volt) for plot, RESUME 
2. enter y scale (in/volt) for plot, RESUME 
3. depress FLAG if printout is desired 

4. enter X scale (units/5 volts) for printout, RESUME 
5 


enter Y scale (units/5 volts) fcr printout, RESUME 










Sense 
switch 
up? 


depress RESUME 





NO (line plot) 


HALT 
depressed? 






NO 


1 next point is plotted 
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Statement Number 


Octal 


NO fF © 


ON OM A 


1 


NNW Uh WN F&F 


20 


mm WwW Wh 


NO WN 


30 


JI HU B&B W NH HY 


Decimal 


0000 


001 


002 


003 


-~— OO DNA MF UD $$ OW WAA WM FWA ~~ OM @ YX Av & WD - 


APPENDIX D 


DATA REDUCTION PROGRAM LISTING 


Code 


266 
160 
220 
016 
266 
161 
220 
016 
266 
162 
220 
016 
266 
163 
160 
314 
260 
025 
314 
220 
034 
302 
247 
025 
161 
214 
241 
051 
360 
326 
315 
273 


Explanation (if needed) 


load IX register | 

with 160 to transfer 
Jump to octal data 
address 16 | te, pikotter 
load IX with 

iL 251 to transfer 
JUNpPREO! Ce tall Y data 
address = 

load IX with 

hee pen up 
jump to octal 

address 

load IX with pen down 
ies 


store IX into IRO 
OPCB (output control byte) 
if 2h of IX # O jump 

to octal address 25 

OPCE 

Jump to octal 

address 34 

RINS (read input status) 
if Hotere ce— 1, 
octal address 25 


jJunp to 


EecallPirOtreqister to 1x 
OPCB 
iz L2 of IX = 1, 


octal ad@ress 51 


jump to 


clear IX 

transtemedigit of E to IL 
OPDB (sign transfer) 
load@oiigngeh Lil Lil 
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40 


“OO Mm &B WN F 


5 


© 


NO Mm & W DN F&F 


6 


oO 


“WO WM B&B WD F 


70 


“NO OW B® WH -F 


100 


“a Om & W DH F 


© 
© 
iol 
Le) 


004 


005 


006 


007 


~OWO DN AM FWN | OW WNAA UM FWA —" OU WA AU FWA OW @ORQIAW & uw 


277 
322 
315 
336 
322 
315 
336 
$22 
315 
057 
066 
052 
114 
000 
063 
377 
240 
120 
241 
ag 
242 
135 
243 
143 
244 
151 
245 
156 
246 
163 
247 
170 
111 
000 
024 
002 
005 
006 
020 
057 
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load DL wath O11 111 111 


- load IX with (DC)* 


OPDB (tramsfer exponents) 

decrement DC by 1 step 

load IX with (DC) * 

OPDB (tramsfer first 2 digits) 

decrement DC by 1 step 

load IX with (DC)* 

OPDB (tramsfer next 2 digits) 
RESUME (emd plot subroutine) 
sylibolic 52: subroutine to 

| ee binary to decimal 
clear data register 

O 

clear entry register 

no operation 

jump if L1 = 1 to 

octal addres 120 

Jump 12 ze = 1 to 

octal address 127 

jump if L@4 = 1 to 

octal address 135 

jump if L& = 1 to 

octal address 143. 

jump if H&F= 1 to 

octal address 151 

junp if H2 = 1 to 

octal address 155 

jump if H4 = 1 to 

octal address 163 

jump if H&® = 1 to 

octal address 170 


to recomvert 
input datta 

E register now 
contains input 
data valwme 


RESUME 


Teh, 


; Pet Titi 








110 


~ 
NO 
NF ON WD MO BW Db & 


! 
mW & WwW 


~] 


130 


fb 
q> 
SO OW B® WDNH FON DM BW ND FF 


jan 
WN 
© 


“a HWM B&B W NH FH 


oO 
Oo 
~ 
Le) 


008 


009 


010 


O11 


~SCODANANF UND HM OU MAAN FWA HK OWUDRAAWNF® WA — 


OwoON A VI & UI 


377 
BG 
377 
377 
377 
377 
377 
377 
001 
002 
010 
113 
000 
220 
062 
006 
004 
113 
000 
220 
064 
003 


002 


113 
000 
220 
066 
001 
006 
113 
000 
220 
070 
010 
113 
000 
220 
07.2 
004 
113 


no operations 


4 dd 128 to main data register (MDR) 


2 O if most significant 
8 but (LL) O82 x recister 
+ is on (Ll = 1) 
Q 
jump to octal 
address 62 
6 add 64 to MDR 0 
é Ie) ee 
+> 
Q 
jump to octal 
address 64 
J 
2 add 32 to MDR O 
+ if L4 =1 
0 
jump to octal 
address 66 
VA 
6 add 16 “toe HDR 0 
+ if LS = 1 
4) 


jump to octal 


8 address 70 
+ add 8 to iiDR 0 
jump to octal 
R° address 72 


4 
+ Wen: EOvMDR TO 
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160 


“NO WB WH fH 


17 


© 


“SO WO ® WN & 


200 


“NO UO BP WN & 


210 


N 
Wy) 
“NO uw BR WH FON DU BW NE 


© 
=e. » 


012 


013 


© 
> 


015 


WD 


—-~OWODONA MW EFUND | OW DAAUMN FWD KH OW DWNA WM FWDND = OW @D AAW & w 


000 
2240 
O74 
002 
113 
000 
220 
Gitis 
001 
113 
000 
220 
100 
377 
377 
377 
377 
266 
000 
314 
305 
163 
266 
000 
314 
305 
057 
066 
075 
200 
014 
200 
010 
056 
057 
066 
070 
023 
111 
002 


if H2 = 1 
jump to octal 
address 74 
add 2 to MDR O 
if H4=1 


jump to octal 
address 76 
add 1 to MDR O 
iz H8 = 1 


jump to octal 
address 100 


no operations 


next is channel select and 


data transfer code ror A/D 
converter 


load IX with 

(code for Y data channel) 
OREGB) get y value in IX 
IPDB) 

store. IX in IR] register 
load IX with 

(code for x data channel) 
OPCB)get x value in IX 
IPDB) 

RESUME 


pre: Wis) elasyghe 
2 


plot subroutine 
branch te soca: 
address 14 
branch to.oOta:| 
address 10 

hale 


RESUME 


L Fiat Sin x data scaling 


and transfer subroutine 


EX MDR 2 for x scale 





TU UBL 





Doolir ua 


Bie on 





230 


24 


© 


“a OO OM Rh WD 


25 


© 


“OO OM Rf WD EF 


260 


“a OF Om B&W DY FF 


27 


© 


~ DW OO B&B WN 


“SN OO MW & WH 


© 
uw 


016 


017 


018 


019 


RO 


—~O WOO NRMA NSFW ND |" OWAANA MN FWHM |" OW DNA MN FUND |= OW WHOA WN PP Ul 


021 
111 
010 
020 
220 
000 
066 
071 
023 
111 
005 


(021 


111 
O11 
020 
220 


004: 


066 


005 


065 
166 
117 
005 
200 
010 
056 
110 
002 
060 
056 
110 
003 
060 
056 
110 
004 
060 
056 
110 
005 


add x origin 


4 


jump to x data entry of 
O plot routine 


- arc sine Y data scaling 


and transfer subroutine 


E X MDR 3 for Y scale 


add Y oriqin 


advance paper 


clear flag 


} jump to Y data entry 


4 plot routine 
symbolic 5 scale factors 
Ss storing subroutine 


set decimal point to 


J 3 places 


Jump to pen up entry 


y | O= DICE eoutine 
halt 


2 enter and print 
£Or pilot 


hale 


enter and print 
JS for printout 


halt 
enter and print 
y EOL splot 
halt 
fi eich pis 3 mit 
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x 


x 


VW 


Vv 


scale 


scale 


Y scale 


scale 





300 


“NO OW Rh WD & 


310 


“NI OM B&B WN FPO NW MO BW NFP ON DUM RW YD & 


340 


“NO MW & W DN 


060 
065 
066 
006 
200 
201 
200 
054 
127 
016 
067 
074 
127 
067 
070 
165 
200 
054 
127 
016 
067 
054 
127 
067 
071 
127 
023 
067 
O75 
126 
067 
006 
066 
074 
110 
001 
11] 
023 
004 
266 


ror printout 
advance paper 
symbolic 6 
S main routine 
branch to octal address 201 
yan x and y data 
branch to octal address 54 
Ye } 


evaluate x data 


branch to symbolic sin 
(x scale and plot routine) 


by branch i£ flag is set to 
symbolic r(x print 
routine 
r 
Br 
JU 


recall IRL register to IxX(Y data) 


% fos magdata 

branch of flat is set to 
symbolic 1/x (Y print 
routine) 


branch to v scale and 


at 


plot routine 


jump to swmbolic 6 again 
(start of main proaran) 


S | 
br 
X 
Z 
] 
6 
symbolic xr 
hs print routine 


, save x value 


scale x 
load IX with 
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350 


N HAM hh WwW ND F 


36 


© 


NO OBR WD 


37 


> 
© 
“SO MW BPW DNO FON BD MN hm WW NEF OC 


> 
~ 
© 


“~ HO WM BW DN FF 


© 
Le) 
Li 
RO 


O24 


O25 


026 


O27 


—~—OCWODANWM FUND KH OW DADAM FUND HH OCW MAAN FUDYD—OWDRAAWe WY 


354 
150 
111 
001 
OE a 
066 
054 
110 
004 
111 
023 
005 
266 
3515 
150 
176 
111 
001 
O57 
066 


012. 


062 
060 
200 
000 
200 
004 
056 
001 
000 
060 
200 
000 
200 
004 
O66 
000 
056 
110 
010 


0 
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symbolic 
and scale check routine 


11 101 100 
print scaled x and x 


return x value 
RESUME 


Bare 1/x 


Y print routine 


same Y value 


scale Y 
load IX with 
LL Oi rer 


print scaled Y and Y 
print data 


return Y value 
RESUME 


* plotter "“zerc 


send 0 to plotter 


send 10 to plotter 


grid coordinate routine 


} symbolic 0 


enter X grid coordinate 





© 
RO 
~) 
8) 


420 060 
200 

000 O 

056 enter Y grid coordinate 
110 ¢ ; 
O11 9 

060 

200 

004 4 

066 } symbolic 1 axis plot 
00] y, routine 

056 enter x iiaximum 


110 ‘ 

000 0 

060 

056 enter x minimum 
110 ‘ 

001 , 

060 

056 enter y maximum 
060 

110 ’ 

004 - 

056 enter y minimum 
060 

110 é 

005 s 

066 symbolic 2 
002 2 

056 enter x scale 
060 

110 } 

002 > 

056 enter y scale 
060 

. 110 ’ 

003 J 

111 ' 

000 0 

111 ¢ 


OWUODRAN WM ew 


028 


Aas 
W 
SND OM & WN FON DM & WH 


440 


029 


“SNH OM & WD 


45 


030 


“ND MW B&W ND FH 


460 


O31 


© 
~S9UODAAMFWNH KH OW MAAN FWHD | CWO W@W ADAM eB wD = 


“SO OM & W DN FE 
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470 


: Sa) 
© 
NO MW hWDNYD FP ON DOD UO hW DY 


510 


“NO Ul BR WN F&F 


52 


Oo 


NO OB WN 


53 


© 


“4 OO BR W NY FF 


© 
Yad 


O52 


033 


O34 


O35 


—~OWODNAWMF UND HS OWODAAM FUND HS OWUDWAAWM FWA HOU D AAW &w 


A) 


022 
001 
353 
111 
024 
002 
060 
111 
004 
111 
022 
005 
353 
111 
024 
003 
060 
056 
060 
110 
006 
056 
060 
110 
007 
066 
003 
062 
111 
001 
O41 
127 
067 
070 
000 
127 
067 
O71 
127 
067 


t 
+ 
@ 
compute and print length 
f of xX axis 
4 
¢ 
Ss 
t 
7 
J 
compute and print length 
of y axis 
enter interval between graduations 
4 Ortx 
c 
enter interval between 
graduations o£ y 
; 
7 


symbolic 3 | 
J axis plot routine 


¢ 
7 
€ 
S transfer x value 
0 
t transfer y value 
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540 


nn 
Ov 
NO MO BWNYN FP ON DMN B&B WN rF ON DBD MOB WN FY ONY DO WM Rh WD F& 


600 


VN HO MW RW YN FY 


0352 


036 


037 


038 


O39 


OM oO NA WY & w 


—~OWODNAMNFUN KH OW MANAUS UNHD HOW DAAWNWeE WA — 


077 
066 
050 
043 
127 
067 
070 
000 
127 
067 
071 
043 
127 
067 
070 
012 
001 
111 
021 
O11 
200 
004 
043 
127 
067 
070 
012 
001 
111 
022 
O11 
013 
200 
004 
043 
Wea 7 
067 
070 
000 
127 


t symbolic lg 


routine to plot x axis 


oO 





eipe ne 








610 


NO MO” ff WN F 


62 


© 


&m Ww NF ON DO MO Fh WwW NY EF 


“a HO MO BPW NH FON BD OB WN FON HD WN 


© 
Yl 
"@) 


040 


O41 


042 


043 


—) 


—~AOODNAWMF WD HK OW DWAAMFWDHD | OW WA AY & WD OO WD AA YW & Ww 


067 
071 
111 
006 
041 
111 
000 
022 
043 
020 
126 
021 
067 
050 
126 
020 
067 
050 
200 
010 
000 
127 
067 
070 
062 
111 
005 
04} 
127 
067 
07 1 
127 
O67 
O77 
066 
051 
000 
127 
067 
070 


ay 


by 


¢+f& usd 3 + 


ue 


~) 


symbolic e” 


y axis plot routine 


on 








—_ 





: 






‘ oF —» 5 ‘ Ss» .s @ 5 “a 


imately iiiiel!- ot 
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TABLE VII 


DATA REDUCTION PROGRAM SCRATCH PAD REGISTER USAGE 
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